Solitary waves and yrast states in Bose-Einstein condensed gases of atoms 
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Considering a Bose-Einstein condensed gas confined in one dimension with periodic boundary 
conditions, we demonstrate that, very generally, solitary- wave and rotational excitations coincide. 
This exact equivalence allows us to establish connections between a number of effects that are present 
in these two problems, many of which have been studied using the mean-field approximation. 
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Introduction. Remarkable experimental developments 
in the field of cold atomic gases now permit the realiza- 
tion of systems not previously accessible. One example 
is the recent fabrication of toroidal traps that realize pe- 
riodic boundary conditions [ll-dl- Tight, elongated traps 
now permit the study of quasi-one-dimensional motion 
since the energy associated with transverse excitations is 
much higher than the interaction energy ^ . A very tight 
toroidal trap is expected soon and will permit the study 
of periodic motion in a quasi-one-dimensional system. 

Such advances will make it possible to study many of 
the interesting phenomena predicted by the various non- 
linear models of such systems that have been studied in 
recent decades. In the case of bosonic atoms, for exam- 
ple, even the ground state is nontrivial. When the ef- 
fective interaction is attractive, a localized density max- 
imum forms for sufficiently strong coupling 0, Q • When 
the effective interaction is weakly repulsive, the density 
is homogeneous. As the interaction strength increases, 
however, the gas eventually enters the Tonks- Girardeau 
state of hard-core bosons, in which the bosons resemble 
fermions in many respects [l^ . 

Another interesting aspect of bosonic atoms confined 
by a ring potential is their excitation spectrum pl] - [l3| . 
Two fundamental forms of excitations, solitary-wave ex- 
citations and rotational excitations, have been investi- 
gated. In the first case, one is interested in traveling- 
wave solutions for which the wave propagates around the 
ring with a constant velocity and without change of form. 
In the second case, one is interested in the lowest-energy 
state of the system for a given angular momentum, which 
is the so-called "yrast" state. 

Clearly, both are excited states of the system, and it 
is interesting to investigate if and how they are related. 
Bloch's theorem [3] implies that states of different val- 
ues of the winding number are connected via excitation 
of the center-of-mass. This implies that the dispersion 
relation, i.e., the energy of the system as function of the 
angular momentum, is a periodic function on top of a 
parabolic envelope function due to center-of-mass excita- 
tions. The yrast state is (by definition) the lowest-energy 
state for a fixed angular momentum. Therefore, the naive 
expectation is that the dispersion relation appropriate 
for solitonic excitation should have a structure similar 



to that of the yrast spectrum, but with a higher energy. 
However, in a recent study jlS], Kanamoto, Carr, and 
Ueda have provided numerical evidence that "quantum 
solitons in the Lieb-Liniger Hamiltonian are precisely the 
yrast states" . The main result of the present study is to 
demonstrate analytically that solitonic and yrast states 
are identical. 

The present study contributes to the long-standing 
problem of the excitation spectrum of a one-dimensional 
Bose gas with periodic boundary conditions. In his sem- 
inal work Lieb demonstrated that this spectrum has 
two branches. The usual Bogoliubov spectrum which is 
linear for long wavelengths and quadratic (i.e., single- 
particle like) for higher energies. It is shown here that the 
other branch, which has been identified as corresponding 
to solitary- wave excitations [l^ [l3| , can also be regarded 
as the lowest-energy state of the gas for a fixed value of 
the angular momentum. 

In the following we first demonstrate the equivalence of 
yrast states and soliton states for the case of an axially- 
symmetric ring. This result is completely general and 
makes no assumption about the explicit form of the 
many-body state. We then examine the limit of weak 
interactions, where the mean-field approximation is an 
excellent description of both states, and discuss the links 
between various effects that appear in these two seem- 
ingly different problems. 

The equivalence of rotational and solitonic excitations. 
Let us start with the yrast states. For a given Hamilto- 
nian, if, the constraints of a fixed particle number and 
a fixed angular momentum can be imposed with the in- 
troduction of two Lagrange multipliers, /i and fi, which 
can be interpreted physically as the chemical potential 
and the angular velocity of the trapping potential. (We 
assume here that the gas has equilibrated in the rotat- 
ing trap). If ^j,(a;i, X2, ■ . ■ ^xn) is the yrast many-body 
state, variations of the energy functional. 



E{<i>y,%) = j %H^ydXl...dXN 

/ -^*y-^ydxi...dxN I '^*yL^ydxi...dxN, (1) 
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with respect to 4'* yield 



- - f^i*,y = 0. 



(2) 



Turning to the solitonic solutions, we assume that 
these have the form 

^sixi,X2,..., XN,t) = *s(2i,-Z2,...,-Z7v)e-'''*/'', (3) 

where Zi = Xi — ut, where u is the velocity of propagation. 
Since iTid^s/dt = H<^,, 



H^>s - M*s - uP^s = 0, 



(4) 



where P is the momentum operator. Comparison of 
Eqs. © and ^ reveals that and satisfy the 
same differential equation and the same boundary condi- 
tions. They are thus identical with the obvious equality 
u = riR, where R is the radius of the ring. 

This equivalence is extremely general and requires only 
that the Hamiltonian is invariant under the transforma- 
tion Xi Zi. While the kinetic energy and the two-body 
interaction evidently meet this requirement, this is not 
always the case in the presence of an external one-body 
potential, V{x). The above results apply when V{x) is 
axially symmetric and time-independent. The imposition 
of strict axial symmetry has important effects on the ex- 
act wave functions and It forces them to have 
an axially-symmetric single-particle density distribution, 
and it implies that the condensate is, in general, frag- 
mented. However, these conclusions do not necessarily 
reflect the behaviour of real physical systems. In practice, 
it is impossible to avoid weak anisotropics in the trapping 
potential. Such symmetry-breaking anisotropics, even 
those which are vanishingly small in the limit of large 
atom numbers, are sufficient to break the axial symme- 
try of the single-particle density and to restore an un- 
fragmented condensate [20h23i . 

In the following, we focus on the mean-field approxima- 
tion, which is valid for sufficiently weak interactions. At 
mean-field level, we seek yrast states with a constrained 
average value of L, rather than insisting on having eigen- 
states of L. It should be noted that although the yrast 
state and its corresponding soliton state are rotationally 
symmetric, this rotational symmetry breaks within the 
mean-field approximation for both of them. It has been 
shown, however, in the limit of weak interactions that the 
energies of the axially-symmetric and broken-symmetry 
yrast states are identical to leading order in with dif- 
ferences of order 1/N [H [H. In the case of solitonic 
solutions we seek travelling-wave solutions, which prop- 
agate with a constant velocity. As shown more gener- 
ally above, both problems reduce to the same differential 
equation for the order parameters, with the result that 
i^y = i^s = V'- Specifically, 



ifiu—— = -—r 

dx 2M dx^ 



(C/o|^|'-m)V. 



(5) 



In this equation we have assumed contact interactions, 
where Uq is the matrix element for zero-energy elastic 



atom- atom collisions (in one dimension), and where u = 
ilR as above. We now examine the consequences of this 
equivalence. 

Consequences of the equivalence. The solitonic solu- 
tions for a Bose-Einstein condensate confined to a ring 
potential with a repulsive effective interatomic interac- 
tion, Uq > 0, are known to be Jacobi elliptic functions 
[H, I23 with density 



n{z) 



I sn 



K{m)z 
ttR ' 



(6) 



where sn(a;|m) are the Jacobi elliptic functions, rimin and 
ri,„ax are the minimum and maximum values of the den- 
sity, and K{m) is the elliptic integral of the first kind. 
The parameter m is determined by the equation 
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where ^0 is the coherence length corresponding to a den- 
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The velocity of propagation is given by 
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Here c is the speed of sound in a homogeneous gas of 
density nmax, Afc^ = f^maxf^O: and q is the phase winding 
number. According to the above argument, the yrast 
states in mean-field approximation are given exactly by 
the same expressions. 

Dispersion. We turn now to the dispersion relation, 
i.e., the energy as function of momentum, appropriate for 
both solitons and yrast states. Given the solitary-wave 
order parameter in the form ips = y^e^'^, one can deter- 
mine the energy and momentum per particle according 
to the following formulas: 



N 



and 



ttR 



-itR 




(9) 



p^{h/N) I n{d(i)/dz)dz. 

-ttR 



(10) 



The analytic form of 4>{z), which is somewhat compli- 
cated, is given in Ref. [231. The desired E{p) follows im- 
mediately given E{u) and p{u). Figure 1 shows the re- 
sults of this calculation for two different situations. Only 
the first two branches of this function are shown. Note 
that E{—p) = E{p). Clearly, the dispersion relation has 
the structure expected from Bloch's theorem as stated 
above. 

It is convenient to introduce the dimensionless quan- 
tity 7 = NUqM R/ (irh^), which is the ratio between 
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the interaction energy noUo, with no being the homo- 
geneous density rig = N/{2ttR), and the kinetic energy, 
Eq = /{2MR^). In Fig. 1 the value of 7 has been 
chosen to be 120/7r^ in the upper curve, and 3/(47r^) 
in the lower curve. In this plot the energy is mea- 
sured in units of Eq = fi^ /{2MB?'), and the momen- 
tum is measured in units of po = R. When p = 0, 
the energy E{p = 0) conies purely from the interac- 
tion energy, which is UqN / {AnR) or 7/2 in units of 
Eq. The maximum value of the momentum per atom 
in each "quasi-periodic" interval is 27rpo- Finally, the 
energy difference E{p = 27rpo) — E{p = 0) = Eq, and 
E{p = 47rpo) - E{p = 2'kpq) = (2^ - l)£'o. More gen- 
erally, these energy differences are all the odd multiples 
of Eq, due to the excitation of the center of mass, as 
expected from Bloch's theorem. 

In the physics of solitary waves, it is well known that 
the velocity of propagation, u, of a solitary wave sat- 
isfies the equation u — dE/dp. This equality can also 
be viewed from the point of view of the yrast states in 
mean-field approximation. Assume that E{1) is the yrast 
energy per particle as function of the angular momentum 
per particle I and that E' is the energy in a rotating frame 
of reference for which the density distribution of the yrast 
state is stationary. Given that E' = E — lil and that 
dE'/dl = 0, it follows that fl = dE/dl or u = dE/dp, 
which is the formula given above. 

In Ref. (27} it was argued that, when the winding num- 
ber q is zero, the velocity of propagation of a solitary 
wave cannot vanish. This is obvious from the dispersion 
relations of Fig. 1. The region from p = {) to p — irpQ 
(given by the left solid points in the two plots of Fig. 1) 
corresponds to q = 0. The smallest possible value of u 
(for q — 0) is attained at p = irpQ where its value is 
u — fi/{2MR). At this value of p the density develops 
a node, and a "dark" solitary wave forms. The density 
is relatively insensitive to changes in the momentum in 
the vicinity of this point, but there is a violent and dis- 
continuous change in the phase of the order parameter 
that is associated with the change of the winding num- 
ber. More generally, the solitary waves become dark and 
the winding number changes whenever p is an odd mul- 
tiple of TipQ. The first two such points, where q changes 
from to —1 and from —1 to —2, are given as solid points 
in the two plots of Fig. 1. 

It is also interesting to note that the interaction energy 
(which scales like N/R) decreases relative to the zero- 
point energy [which scales like Tt? /{MR^)] as the radius 
of the ring decreases for a fixed value of N/R. The ratio 
of these energies is the parameter 7, which is cx [N / R)R^ . 
As R decreases for fixed N/ R, the interaction energy be- 
comes less important, and the gas approaches the limit 
of a non-interacting system. In the limit of small R, (i.e., 
7^1), the quasi-periodic intervals of the dispersion re- 
lation become linear with a discontinuous first derivative 
between them. The velocity of propagation of the soli- 
tary wave thus approaches constant values, given as odd 
multiples of m = h/{2MR), set by the non- interacting 
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FIG. 1: The dispersion relation E{p), for 7 = 120/7r^ (upper 
plot), and 7 = 3/(47r^) (lower plot). The energy is measured 
in units of Eq = /{2MR^), and the momentum in units 
of po = Ti/R. In the two curves both N/R and Uo have the 
same value, however in the lower curve the radius is 40 times 
smaller than in the higher one. 



problem [27| . This effect is clearly seen in the lower plot 
of Fig. 1. 

From the above remarks it is obvious that, as R de- 
creases, there is a minimum value of the radius below 
which u cannot vanish (26| . The condition for this is 
= ■^6/2, where ^ is the coherence length correspond- 
ing to a density uq = N / {2ttR) . In the language of yrast 
states, this says that there is a critical coupling of 7 = 3/2 
for having a metastable minimum at L/N = I = Ti (or 
equivalently at p = 2iTpo [3, Q. For smaller values of 7, 
the slope of the dispersion relation is positive as Z — ^ 7i~ 
[or p (27rpo) ], and the possible metastable minimum 
at this point disappears. 

Energetic stability of the solitary waves. The creation 
of solitary waves in ring-like potentials will be of consid- 
erable interest. The intimate connection between soli- 
ton and yrast states discussed here enables us to draw 
a number of conclusions about the possible generation 
of solitary waves and their energetic stability. The most 
obvious of these is that the lowest-energy state with a 
given (L) will also be a solitary-wave state. Ideally, one 
would like to observe a family of solitary-wave solutions 
much like those observed in harmonic traps [28, 29], i.e., 
density depressions travelling with a velocity which is less 
than that of sound and which decreases as the depression 
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becomes deeper. 

In the case of a repulsive effective interatomic interac- 
tion, the curvature of the dispersion relation is negative 
(as shown also in Fig. 1) 30]. The "kink" states with a 
momentum per particle which is equal to an integer mul- 
tiple of 27rpo SLie the only possible local minima in the 
absence of rotation of the trap. The occupation of these 
states will give rise to persistent currents. Since these 
states have a homogeneous density distribution and a 
non-zero circulation, they would be difficult to observe 
experimentally [Slj. The same conclusion applies to a 
rotating trap. One way to overcome this difhculty would 
be to consider systems with attractive effective interac- 
tions 15]. In this case the curvature of the dispersion 
relation is positive 's^, and persistent currents cannot 
be stable. In contrast, solitary waves can be stabilized 
by rotating the trap. The energy in a frame rotating with 
angular velocity is E' = E — in, and E' can have local 
minima for a continuous range of I which is determined 

by n. 

Summary. In the present study we have presented a 
unified picture of two seemingly different forms of exci- 



tation of a Bose gas moving in one dimension with pe- 
riodic boundary conditions, namely rotational and soli- 
tonic excitation. We have shown that yrast states and 
solitonic states are identical whenever the Hamiltonian 
is axially-symmetric. This result is a simple consequence 
of the more general observation that the problem of find- 
ing wave functions that propagate with a fixed velocity 
is equivalent to the problem of minimizing the energy 
subject to a constraint of fixed angular momentum. The 
only necessary condition for this equivalence, which holds 
for all functional forms of the wave function, is the ro- 
tational invariance of the Hamiltonian. This equivalence 
permits the consideration of such systems from two dis- 
tinct points of view, and we have noted certain aspects 
of a unified picture that result from this observation in 
the mean-field regime, which is valid for sufficiently weak 
interactions. 

This connection between solitons and yrast states has 
been suggested previously [15] and the arguments offered 
here are admittedly elementary. It would appear, how- 
ever, that the rigorous equivalence demonstrated here is 
neither generally known, nor adequately appreciated. 
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